function [Energy, DOS] = ExtendDOS(E, D)
he = E(2)-E(1);
Ne = length(E);
padNum = round(Ne/2);

% Energy Grid Extension
Energy = linspace(E(1)-padNum*he, E(end)+padNum*he, Ne+2*padNum);
% DOS Extension
DOS = padarray(D, [0,padNum]);


end